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We have measured the hyperfine structure and isotope shifts of the 402.0 nm and 399.6 nm res¬ 
onance lines in ^^®Th"''. These transitions could provide pathways towards the excitation of the 
229rph low-energy isomeric nuclear state. An unexpected negative isotope shift relative to ^32^^]^+ 
is observed for the 399.6 nm line, indicating a strong Coulomb coupling of the excited state to the 
nucleus. We have developed a new all-order approach to the isotope shift calculations that is gener¬ 
ally applicable to heavy atoms and ions with several valence electrons. The theoretical calculations 
provide an explanation for the negative isotope shift of the 399.6 nm transition and yield a cor¬ 
rected classification of the excited state. The calculated isotope shifts are in good agreement with 
experimental values. 


The existence of a low-energy nuclear isomeric tran¬ 
sition in 229Th [Tf[5] has stimulated the development 
of novel ideas and concepts in the borderland between 
atomic and nuclear physics HIS], with conceivable ap¬ 
plications like nuclear clocks [SHin] and 7 -ray lasers 
m- The wavelength of this very narrow nuclear tran¬ 
sition is predicted to be 160(10) nm [3]. The big 
uncertainty in the transition wavelength, its location 
in the vacuum ultraviolet range, and the radioactiv¬ 
ity of 229Th with a halflife of about 7900 years present 
experimental difficulties that have so far prevented op¬ 
tical spectroscopy of the transition. A promising ap¬ 
proach towards an excitation of the isomer is the use 
of electronic bridge or NEET processes (nuclear excita¬ 
tion by electron transition) |12lHl] . that involves find¬ 
ing sufficiently strong electronic transitions close to the 
nuclear resonance. The Th^ ion, with its large den¬ 
sity of states, offers a high probability of locating such 
electronic bridge pathways. Using the long-lived 
isotope we have investigated two-photon excitation of 
electronic levels of Th+ in the expected energy range 
of the isomer and have found 43 previously unknown 
levels m- 

A well known influence of the nucleus on electronic 
spectra of atoms is the isotope shift (IS) that for heavy 
atoms is dominated by differences in the overlap of nu¬ 
clear and electronic charge distributions, the so called 
field shift (FS). Extensive studies of isotope shifts in 
the dense level schemes of Th and Th+ have been per¬ 
formed HMH providing shifts of more than 350 ex¬ 
cited states of 2307’]^-!- respect to 232xh+ [inj. In 
all but 3 cases, the shift of the excited state is smaller 


than that of the 6cP7s ground state (GS), indicating 
the smaller overlap of the excited electronic configura¬ 
tion with the nucleus for almost all states. Although 
a nearly complete compilation of the 15 lowest elec¬ 
tron configurations is published m, only very few low- 
energy excited states belonging to fs^ and s^p config¬ 
urations with strong Goulomb interaction with the nu¬ 
cleus are known in Th+. Such states will be of relevance 
for the implementations of electronic bridge or NEET 
processes. 

Based on laser spectroscopy of trapped ions, we 
present here measurements of the isotope shifts and 
hyperfine structure (HFS) of 229xh~'' for two of 
the strongest Th+ resonance lines at 399.6 nm and 
402.0 nm. A positive IS of the 402 nm line of 229Th''' 
relative to 2327 ^]^+ Qf ^9 GHz has been reported al¬ 
ready [TB]. Unexpectedly, we observe a negative iso¬ 
tope shift (i.e. a lower transition frequency) of the 
O 3/2 25027i/ 2 line at 399.6 nm with respect to 

232Th+, indicating a strong field shift of the excited 
state (states are labeled by their energy in cm“^ and 
the angular momentum J). Since the 25027i/2 state 
was previously identified as 5f6(P jlO], a negative iso¬ 
tope shift would contradict the available IS data for 
Th+ as discussed above. 

Resolving this issue requires precision theoretical IS 
calculations. This is a difficult endeavor due to several 
factors: (1) Th+ is a heavy trivalent system with elec¬ 
tronic configurations containing 5/ electrons that are 
very strongly correlated with core electrons; (2) there 
is a very large number of energy levels between the 
ground state and the 25027i/2 level (it is the 82nd odd 
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level); (3) the state of interest is very strongly mixed 
with nearby states complicating the calculations and 
identification of the levels. Therefore, a very accurate 
new approach is needed for the Th+ IS calculation. 

We have developed a method to calculate IS for 
multivalent heavy systems based on the hybrid ap¬ 
proach that combines configuration interaction and 
all-order linearized coupled cluster methods (Cl+all- 
order). While the Cl+all-order method has been ap¬ 
plied to calculating energies, transitions properties, and 
polarizabilities |21Il23] . this is its first implementation 
for IS calculations. This is also a first application of the 
Cl+all-order method for such highly-excited states of 
multi-valent atoms for any system. An ability to accu¬ 
rately calculate IS for heavy ions opens a pathway for 
the accurate determination of nuclear radius differences 
between isotopes of heavy elements based on ab initio 
calculations. 

In our experiments we use a linear Paul trap, de¬ 
scribed in Refs. [ISIEI]. The trap is loaded with > 10® 
Th+ ions by laser ablation from a dried Th(N 03)4 solu¬ 
tion deposited on a tungsten substrate. Collisions with 
argon buffer gas at 0.1 Pa pressure cool the ions to room 
temperature and quench the population of metastable 
states that are optically pumped by the laser excitation. 
Pulsed excitation of the 248745/2 and 25027i/2 states is 
provided by a CW extended-cavity diode laser (ECDL) 
with an output power of 25 mW passed through a fast 
acousto-optical modulator. The excitation pulse du¬ 
ration is 70 ns and the pulse repetition rate of 1 kHz 
is adapted to the collisional quenching rate of about 
«1500 /s pi]. Fluorescence of the excited Th^ ions is 
detected using a photomultiplier (PMT). A fast gated 
integrator is used to evaluate the PMT signal during a 
detection window with a duration of about 75 ns. The 
signal is then integrated over several hundred pulses. 
A small fraction of the CW radiation of the excita¬ 
tion ECDL is coupled to a Fabry - Perot interferometer 
used as a reference to measure frequency intervals. Af¬ 
ter loading of Th+ we observe simultaneously the fluo¬ 
rescence signal from both ^^®Th and ^®^Th, indicating 
approximately equal abundances of both isotopes on 
the substrate. When it is necessary to avoid the in¬ 
fluence of ^®^Th^ ions we apply an isotope selective 
resonantly enhanced three-photon ionization to Th^'*' 
m and eject these ions from the trap. 

Determining the IS for ^^®Th with a nuclear spin 
/ = 5/2 requires knowledge of the HFS centroid fre¬ 
quency. The HFS of the O 3/2 —>■ 248745/2 transition of 
229Th+ at 402.0 nm consists of 12 lines with 4 compo¬ 
nents in each of the AF = —1,0,-1-1 branches, where 
F is the total angular momentum. We investigate the 
HFS by Doppler-free saturation spectroscopy and the 
excitation spectrum of both ^^®Th^ and ^®^Th^ iso¬ 


topes is shown in Fig. Fig. [^3 shows the HFS 
of ^^^Th"*" with high resolution. The width (FWHM) 
of the resonances is >70 MHz and the width of the 
Doppler broadened lines is r:!760MHz. In Fig. & nine 
out of twelve possible resonances are resolved. The res¬ 
onances caused by F = 1 —)■ F' = 2 and 2—3 1 tran¬ 
sitions are overlapping within the width of the Lamb 
dips. The deformation of the Doppler shape caused by 
component 3 —3 2 is also partially observable on the 
graph. The weak HFS component 4 —3 3 is not ob¬ 
served due to the limited signal-to-noise ratio. 

The HFS of the 399.6 nm line O 3/2 —3 25027i/2 is 
shown in Fig. |^. The 25027i/2 state of ^^®Th+ con¬ 
sists of two components with F = 2,3, therefore the 
HFS of the transition comprises 6 lines. The hyperfine 
components of appear at frequencies below the 

position of the ^®^Th“'' line. To exclude a possible mask¬ 
ing of a hyperfine component with the ^®^Th^ Doppler 
broadened line we ejected ^®^Th~'’ with the technique 
mentioned above. This experiment shows that there 
is no hyperfine component of the ^^^Th'*' overlapping 
with the 232rp]^-i- saturation resonance of the 

transition 2—3 3 has a very weak intensity and is not 
observed in the experiment. The resonances caused by 
transitions 3—3 3 and 3—3 2 are resolved by scanning 
with longer integration time (see the inset in Fig. 2b). 
The dip in the middle correspond to a cross-over res¬ 
onance of V-type with the common ground state level 
F = 3 (see Fig. |^). 

In order to determine the hyperfine splitting factors 
A and B of the 402.0 nm and 399.6 nm lines we measure 
the frequency intervals between the Lamb dips relative 
to the transmission peaks of the interferometer. For 
the 402.0 nm line the measured frequency intervals are 
fitted in a least squares adjustment to determine the A 
and B factors for the ground and the 248745/2 states. 
The algorithm attempts fits for all viable assignment 
combinations of F values for the nine resolved tran¬ 
sitions shown in Fig. The result is shown in Ta¬ 
ble m The uncertainty results from inaccuracy of the 
measured frequency intervals, including a small contri¬ 
bution from the determination of the reference interfer¬ 
ometer FSR. The HFS factors for the ground electronic 
state (GS) in ^^®Th’*’ are already known with high ac¬ 
curacy [IS] (see Table |T| and our results are in good 
agreement. 

For the 25027i/2 state the B factor vanishes because 
of the absence of a quadrupole moment for J = 1/2. To 
determine the A factor of this level we fix the A and B 
values for the GS. We use the same algorithm to fit the 
frequency intervals between resolved transitions of this 
line. The uncertainty of the A value of the 25027i/2 
state is determined by the accuracy of the frequency 
intervals and the uncertainty of the GS factors. Due 
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Laser detuning (GHz) 

FIG. 1. a) Fluorescence signal of the 402.0 nm line of 
229 tjj-i- 232rpj^-i- jgotopes. b) HFS signal of the 402.0 nm 

229Th+ line. The detuning counts from the recalculated 
HFS centroid position of the ^^^Th^ isotope line. The com¬ 
ponents are labeled with the total angular momenta of lower 
and upper hyperfine levels. The dashed arrows in (b) show 
the recalculated positions of the 3 —> 2 and 4 —> 3 reso¬ 
nances. 

to the weaker signal-to-noise ratio and smaller A factor 
of the excited state, the uncertainty of the fit is higher 
than that for the 248745/2 state. Using the obtained 
values for A and B factors, we recalculate the position 
of the hyperfine components which are not observed. 
The calculated splitting of the HFS components 1 —?► 2 
and 2 —>■ 1 of the 402.0 nm line is 18 MHz, explaining 
the unresolved resonances of these lines. As a test of 
consistency, we calculate the complete shape of the sat¬ 
uration HFS spectra according to Ref. pS] , ht it to the 
experimental signal and find that it confirms the results 
presented above. 

TABLE I. Hyperfine splitting factors of ^^^Th^ electronic 
levels. The second row for the ground state shows the val¬ 
ues obtained in Ref. These results represent the fit 

where the A and B parameters for the ground state are fixed 
to the values from Ref. |25j . 


Level cm ^ 

A [MHz] 

B [MHz] 

Reference 

0 

-444.2(3.4) 

-444.2(1.9) 

308(13) 

303(6) 

This work 
Ref. |25l 

24874 

489.2(3.7) 

488.9(3.7) 

-409(18) 

-413(16) 

This work 
This work^“^ 

25027 

-45(19) 

-42(16) 


This work 
This work*“^ 


The IS of 229Th+ is determined from the hyperfine 
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FIG. 2. a) Fluorescence signal of the 399.0 nm line of 
229Th+ and isotopes, b) HFS signal of the 399.6 nm 

229 Th+ line. The detuning counts from the recalculated 
HFS centroid position of the 229rpjj-i- jgotope line, and the 
dashed arrows show the recalculated positions of the reso¬ 
nances 2 —>■ 3 and 3 —>■ 2. The inset shows resolved Doppler- 
free resonances of the components 3 —>■ 3 and 3 —> 2 and 
the V-type cross-over resonance. 


centroid, corresponding to the position of the line for 
A = B = 0. For the 402.0 nm line, the centroid is « 
810 MHz red-detuned from the position of the 3 —>■ 4 
component. The isotope shift of ^^®Th'*’ for this line 
is 18.86(7) GHz, in agreement with earlier, less precise 
measurements [IIlIIH]. The uncertainty is determined 
by the precision of the frequency intervals measure¬ 
ments and the calculated value of the HFS centroid. 
The GS of Th’*’ yields a mixture of electronic configu¬ 
rations with leading components 6cP7s and 6d7s^ while 
the 248745/2 level contains 26% 6d7s7p and 20% 5/6d^ 
|20| . For this transition the FS of the ground state is 
dominant m and determines the positive IS of the 
402.0 nm line. 

The 399.6 nm line shows a negative IS of 
—3.67(5) GHz. The centroid of the line is shifted at 
« 1.5 GHz to the lower frequency from the resonance 
4 —>• 3 (see Fig. The probability to observe a neg¬ 
ative IS in Th'*' for a line originating from the GS was 
predicted to be very small CHI since the s-electrons 
make the greatest contribution to the IS. For the odd 
parity levels a negative IS is observed only for the tran¬ 
sition from the GS to the 83797/2 level which has the 
leading electron configuration /s^. According to j^HI 
the leading configuration of the 25027i/2 state is 5/6ci^. 
The analysis in Ref. m indicates a weak FS for the 
5f6d^ electron configuration which means that the IS 
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TABLE II. Comparison of the g-factors and excitation en¬ 
ergies (cm“^) of Th"*" odd J — Ijl states calculated using 
CI+all-order method with experimental data compiled in 
Ref. 1201 (see also: |29|L The last column gives finite size 
IS constants Kfs in GHz/fm^. All values are counted from 
the GS. 


State 

Ref. [20] Calc. 

J 

g-factor 
Calc. Exp. 

Energy 
Calc. Exp. 

Kfs 

Calc. 

5/6d7s 

5/6d7s 

1/2 

0.243 

0.255 

12390 

11725 

-82.1 

5/6d7s 

5/6d7s 

1/2 

0.558 

0.523 

14805 

14102 

-96.7 

5/6d7s 

5/6d7s 

1/2 

2.577 

2.565 

16207 

15324 

-65.9 

5/6d7s 

5f6d^ 

1/2 

1.618 

1.080 

18307 

17838 

-120.4 

5f6d^ 

5f6d^ 

1/2 

0.405 

1.007 

19309 

18568 

-126.4 

5/6d7s 

5f6d^+ 

5/6d7s 

1/2 

0.730 

0.770 

23456 

22355 

-100.6 

5f6d^ 

7s^7p+ 

6d7s7p 

1/2 

0.568 

0.600 

25445 

25027 

5.8 

5f6d^ 

5f6d^ 

1/2 

2.108 

0.725 

26024 

25266 

-119.7 


of the transition should be determined by the energy 
shift of the ground state which leads to a positive IS. 
A potential explanation for the negative IS would be 
a prior misidentification of the 25027i/2 level and the 
presence of a large admixture of a Is^nl configuration. 

The field shift for isotope A' with respect to isotope 
A is defined as: = ATps S{r'^)^ where iCps is 

the field shift constant and S{r'^)^ = (r^)^ — 

The FS operator, H-pst is a one-particle operator which 
modifies the Coulomb potential inside the nucleus. We 
use the “finite field” method, which means that a per¬ 
turbation is added to the initial Hamiltonian with the 
arbitrary coefficient X: H —>■ H\ = H + AiLps- We find 
eigenvalues E\ by direct diagonalization of Hx and then 
find the FS coefficient of an atomic level as |271 [28] : 


i^FS 


5 dEx 
6 r^ dX 


( 1 ) 


where the nucleus is modeled as a homogeneously 
charged ball of radius r^v. 

We use the Cl+all-order method, with modified 
Hamiltonian given above, to carry out the field shift 
calculation. All computations are repeated three times: 
with A = 0 and A = ±0.01. Then, Eq. Q is used to 
evaluate the field shift constant. The general descrip¬ 
tion of the CI±all-order method and its application to 
the calculation of energies, transition properties, and 
polarizabilities are given in [2IH23]. We have modified 
the Cl part of the method to be able to successfully 
generate much higher excited states than previously 
possible. This advancement will allow further appli¬ 
cations of the CI±all-order method for highly-excited 
states and systems with dense spectra. Another issue 
of the application of the CI±all-order method to IS is 
the requirement of very high numerical accuracy since 


the correction due to IS is small. Therefore, we im¬ 
proved and thoroughly tested the numerical accuracy 
of the approach such as convergence of various itera¬ 
tion procedures. 

While only ns orbitals have significant overlap with 
the nucleus, with pi /2 having much smaller overlap, 
other states will be affected by the inclusion of the IS 
effects due to the self-consistent nature of the poten¬ 
tial. All of the ns orbitals become slightly less bound, 
leading to the expansion of the total self-consistent elec¬ 
tronic potential. This, in turn leads to changes of en¬ 
ergies of the other states as they become more deeply 
bound in the modified potential. The IS of 6d and 5/ 
electrons will be with the opposite sign of the 7s elec¬ 
tron IS. One-electron estimates give the IS of the 5/ 
electron being about 60% of the 7s IS with an opposite 
sign. 

The results of the CI±all-order calculation of the g- 
factors and excitation energies of Th+ odd J = 1/2 
states are compared with experiment [201 [29] in Ta¬ 
ble |TTj We find good agreement for the energy and g- 
factor of the 25027i/2 level of interest. The last col¬ 
umn gives differential finite size IS constants Kps in 
GHz/fm^ for J = 1/2 levels, relative to the ground 
state. The table clearly demonstrates the opposite sign 
of the isotope shift constant for this level in comparison 
to the other states. Our calculations show Is'^lp with 
a contribution of 36% and 6d7s7p with 34% as the two 
dominant configurations for the 25027i/2 level, which 
explains the observations of the negative IS. Contri¬ 
butions from Qd'^7p are 14% and from 5f6d^ 4% only. 
Using the value = —0.33(5) fm^ [251150] for 

the difference in the rms radii, we get —1.9 GHz for the 
CI±all-order IS of the 399.6 nm line, in good agreement 
with the experimental value of —3.67(5) GHz. The diffi¬ 
culty of this particular calculation is the very large cor¬ 
relation correction for the 5/ electrons and very strong 
configuration mixing for the J = 1/2 states. Our value 
for the 402.0 nm IS is 19.0 GHz, in excellent agreement 
with the experimental results of 18.86(7) GHz. 

We have also carried out calculations using 
CI±many-body-perturbation theory (MBPT) [51] pre¬ 
viously employed for IS calculations [55] . The 
CI±MBPT results for the 399.6 nm and 402.0 nm lines 
are 8.1 GHz and 26.0 GHz, respectively, in very poor 
agreement with experiment. This demonstrates the 
success of the CI±all-order method and the need for 
the inclusion of the higher-order effects. 

In summary, we have demonstrated that the CI±all- 
order method can accurately describe isotope shifts in 
such complicated systems as Th+. In turn, the IS mea¬ 
surements in actinides represent excellent benchmark 
tests for the development of these highly accurate the¬ 
oretical methodologies. The results of this work are 
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important for the selection of efficient transitions for 
the excitation of the ^^®Th isomer, as the IS indicates 
the strength of the Coulomb interaction between elec¬ 
trons and nucleus that is relevant for NEET processes. 
The comparison of isotope and isomer shifts will give 
access to the Coulomb contribution to the nuclear tran¬ 
sition energy, resolving a dispute about largely different 
theoretical predictions of the sensitivity of the nuclear 
transition frequency to the value of the fine structure 
constant [55] , 
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